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Abstract 

The conditions required for a cohesive zone model (CZM) to predict a failure 
load of a cracked structure similar to that obtained by a linear elastic fracture 
mechanics (LEFM) analysis are investigated in this paper. This study clarifies 
why many different phenomenological cohesive laws can produce similar fracture 
predictions. Analytical results for five cohesive zone models are obtained, using 
five different cohesive laws that have the same cohesive work rate (CWR-area 
under the traction-separation curve) but different maximum tractions. The effect 
of the maximum traction on the predicted cohesive zone length and the remote 
applied load at fracture is presented. Similar to the small scale yielding condition 
for an LEFM analysis to be valid, the cohesive zone length also needs to be much 
smaller than the crack length. This is a necessary condition for a CZM to obtain 
a fracture prediction equivalent to an LEFM result. 

Introduction 

The origin of the cohesive zone models (CZMs) can be traced back to the “Dugdale-Barrenblatt 
model” [1, 2], The cohesive zone is considered to be a fracture processing zone ahead of the 
crack tip. For most cohesive laws, the traction-separation curves used to model the material 
within the cohesive zone are phenomenological, and hence, are not directly related to the 
physical process in the damage zone which typically is difficult to determine experimentally. 
Regardless, the CZM approach has been widely accepted as a computationally convenient 
fracture analysis tool. If the CZM approach is used in a finite element analysis, the crack 
initiation, growth, and direction of growth can be automatically determined. Many different 
cohesive laws with variances in maximum traction, maximum separation, and shape have been 
proposed; such as the linear softening cohesive law by Camacho and Ortiz [3], the exponential 
cohesive law by Needleman [4,5] and Xu and Needleman [6], the trapezoidal cohesive law by 
Tvergaardand and Huchinson [7], and the polynomial cohesive law by Tvergaard[8], 
Researchers found these cohesive laws generally produce results that correlated well with 
experimental data such as the failure load and crack growth for cracked structures. For linear 
elastic materials, these predictions are comparable to linear elastic fracture mechanics (LEFM) 
results; however, the conditions under which the cohesive zone modeling approach and the 
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fracture mechanics approach are equivalent have not been systematically investigated. Here, the 
equivalence of the two modeling approaches means that they can predict the same failure load 
for a cracked structure. 


The objective of this paper is to investigate under what conditions the cohesive zone model and 
the LEFM approach are equivalent. Linear softening cohesive laws are used in this study. The 
relationship between the CWR (area below the traction-separation curve of the cohesive law) and 
the J-integral value [9] is discussed. Integral equations, for using CZMs to analyze the fracture of 
an infinite plate with a crack under a remote tensile load, are presented [10]. An iterative 
numerical procedure is implemented as a MATLAB® M-file [11] for solving these equations to 
obtain the length of the cohesive zone ahead of the original crack, the opening displacements 
along the cohesive zone, and the remote applied stress at the moment of crack growth initiation. 
The effect of varying the maximum traction of a cohesive law on the predicted cohesive zone 
length and remote applied stress is investigated. The remote applied stress is used in an LEFM 
formula to determine the energy release rate. By comparing the energy release rate with the 
CWR, it is possible to determine whether the cohesive zone modeling approach is equivalent to 
the LEFM approach, i.e. whether it is able to produce a similar failure load for a cracked 
structure to that produced by the LEFM approach. 

CWR and J-integral value 

A CZM is used to study the fracture of an infinite plate with a crack length of 2a subjected to a 
remote applied tensile stress cr x as shown in Figure la. Note that all the equations and analysis 

results presented in this paper are related to Mode I fracture of linear elastic material. In the 
cohesive zones (narrow bands) ahead of the crack tips, the prospective fracture surfaces are 
assumed to be restrained by a cohesive stress that Dugdale took to be the yield stress of the 
material [2]. The traction (cohesive stress) <r and the separation, S = u\ - a 2 , of the upper and 

lower prospective crack surfaces, where u \ and u~ are the opening displacements in the y- 
direction for the upper and lower surfaces, respectively, are shown in Figure lb. 

A CZM model needs a cohesive law to relate the traction to the separation at the same location 
along the x -axis for describing the constitutive behavior of the cohesive zone. The cohesive 
law used in this paper is a linear softening cohesive law as shown in Figure lc, 


cr = cr c (S c -S) / 8 C 0 <8 <8 C 
= 0 5>S C 


( 1 ) 


where <j c is the maximum traction and 8 C is the maximum separation. The cohesive work rate 
(CWR) <D f , or work of separation per unit area, is the area under the linear softening curve 
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( 2 ) 


At the moment of the crack growth initiation, the cohesive zone opening displacement at the 
crack tip reaches maximum separation, 8(a) = 8 C , and the cohesive zone is fully developed as 

shown in Figure 2 [7], During the subsequent crack growth, the length of the cohesive zone and 
the opening displacement along the cohesive zone length are unchanged. The J-integral [9], 
taken along the boundary of a fully developed cohesive zone, the contour T , can be expressed 
as. 


J = \^Wn l -T i -^)ds (3) 

where W is the strain energy per unit volume, n l is directional cosine between the positive side 
(outward) of the normal N and the x -axis, T t is the i ,h component of the traction perpendicular 
to T in the outward direction, u j is the i th component of the displacement, ds is an arc element of 

r. 


For contour T shown in Figure 3, m, = 7] =0and the J-integral value for a fully developed 
cohesive zone becomes [12] 


p 

[a 

8ul 

du 2 

0 

dx 
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(4) 


where J c is the critical J value immediately before the initiation of crack growth and J cn is a 
J -integral value of an infinitesimal contour, not shown in the figure, around the cohesive zone 
tip. For an isotropic material, J czt is related to the stress intensity factor^,, if it exists, at the 
cohesive zone tip 


J 

J CZt r 

E 


(5) 


Since the J -integral is a path independent integral, the J value obtained from any contour 
integral around the cohesive zone such as T 1 in Figure 3 must also equal to J c . If the cohesive 
zone model is equivalent to an LEFM model, the stress field outside of the cohesive zone must 


3 



be near the same as the K-dominant stress field [13] in the LEFM model. In other words, the path 
independent J c can be related to the fracture toughness K c as 


= ( 6 ) 

and the critical energy release rate defined in LEFM as 

G,=J, . (7) 

For cohesive laws that assume the CWR to be the same as the critical energy release rate, 

« \=G C =J C ( 8 ) 

the stress intensity factor at the cohesive zone tip, K czt in Eq. (5), must vanish because J czt must 
be zero to satisfy Eq. (8). Note that G f and K c are material properties and are experimentally 
determined, so is also considered to be a material property. 


Because J czt = 0 in Eq. (4), then J c obtained by the contour around the cohesive zone is the CWR 
® c corresponding to the area under the traction-separation curve. Cohesive zone models with 
the same CWR but different shapes will produce the same J c . If J c is the only parameter to 

determine fracture, it is expected that the same failure prediction can be obtained with different 
shapes of the cohesive law. This may be the reason why CZMs with different shapes can 
generate similar failure predictions. 

Analytical Solutions for Cohesive Zone Models 

CZM analyses of a cracked infinite plate subjected to a remote tensile stress, as shown in Figure 
la, are presented in this section. Five linear softening cohesive laws shown in Figure 4, with 
different maximum tractions but all having the same CWR, were used in this study. The effects 
of maximum traction on the predicted cohesive zone length and remote applied stress were 
investigated. The predicted remote applied stresses were used in an LEFM formula to compute 
the energy release rate at failure. Good agreement between the energy release rate at failure and 
the CWR determines the equivalence of a CZM and a LEFM analysis. 


The analytical procedure presented here follows the paper authored by Jin and Sun [10]. The 
cohesive zone is treated as an extended part of the crack with the stress intensity factor being 
removed at the tip of the cohesive zone. Hence, 



2 p <r(4)dt; 
4nc a s[\ - 1 2 / C 2 


= 0 


( 9 ) 
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where c is the x -coordinate of the cohesive zone tip and er(£) is the traction at location B, as 
shown in Figure lb. 

The total cohesive zone opening displacement 5 at location xis [14] 


<?(*)= ^ Ji ~~\'o(x,4)c7(4)d4 

F \l c nE Ja 


( 10 ) 


where E is the Young’s modulus, and G(x,£) is an influence function given by 

Vl — x 2 /c 2 +Vl-^ 2 /c 2 


G(x,|) = In 


Vl-X 2 /c 2 -y/l^Tc 2 


For the linear softening model shown in Figure lb, Eq. (10) can be expressed as 


( 11 ) 


4ccr, lx 2 4 r_ ~ cr c (5 c -8(%)) 






( 12 ) 


and Eq. (9) can be expressed as 



(13) 


Substituting cr, into Eq. (12), an integral equation is obtained for determining the cohesive zone 
length, p = c-a, and cohesive zone opening displacements (CZODs) £(x) 


nE8\x)<S c 1 


2ccr„ 


c fG(x,4)8(4)d4+ c jl 


f(4) 


% 2 /c 2 




n C lyll-f/c 2 c{ 


(14) 


where S* (x) = 8(x) / d c 


Equation 14 can be solved with an iterative method shown in Figure 5, in which an initial 
cohesive zone length p = c-a is given first to determine the integral interval, [a,c] , for the 
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integral terms in Eq. (14) and then the CZODs along c — a are computed. If the predicted 
opening displacement at the crack tip c) (<a) is not equal to the maximum separation an 

updated cohesive length is given until the S (a) equals^ . The final value of c-a is the 

cohesive zone length p . The remote applied stress <t w can be obtained with Eq. (13), by using 

the final set of CZODs. The remote applied stress is used to compute the energy release rate 
using the following LEFM formulae. 


G(crJ = ^~ (15) 

E 

If the ratio of the computed energy release rate G(cr x ) to the CWR, , is close to unity, then 
the CZM is considered to be equivalent to an LEFM analysis. 

Results 

Analytical results for five cohesive law models with different maximum tractions as shown in 
Figure 4 are presented in this section. All the results are computed at the initiation of crack 
growth when the cohesive zone is fully developed. The effect of maximum tractions on cohesive 
zone length, remote applied stress, the equivalence of CZM approach and LEFM analysis is 
presented. The effect of maximum traction on cohesive zone length can be found in Figure 6, 
which shows that larger maximum tractions result in shorter cohesive zone lengths, p . Two 

curves represent results obtained with two different crack lengths, one has a / l ch = l and the other 
one has a / l ch = 1 0 . Note that l ch is a characteristic length, and in this paper it is defined based on 
Model B in Figure 4, 


‘ch 



(16) 


For the rest of the paper, the term crack length means the half-crack length a in Figure la. The 
two curves in Figure 6 are converged for maximum tractions greater than 2u c , indicating that 

the cohesive zone length may not depend on the crack length for a cohesive law with a large 
maximum traction. However, the cohesive zone lengths depend on crack length for models with 
a maximum traction less than a c . The cohesive zone length for the short crack a / l ch = 1 is 

longer than that of the long crack a / l ch = 10 . 

Cohesive zone lengths are predicted for all five cohesive law models in Figure 4. The cohesive 
zone length reaches a constant value for all cohesive law models as the crack length increases. 
For illustration puiposes, only results of Models A, B, and D are plotted in Figure 7. Models A, 
B, and D have maximum tractions of l/2cr., cr c , and 2a c , respectively. For the largest 
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maximum traction 2 cr, model, the cohesive zone length reaches a constant value at a crack 
length around a = 21 ch . For the second largest maximum traction cr, model, the cohesive zone 
length reaches a constant value at a crack length around a = 5 l ch . The cohesive zone length 
almost reaches a constant value at a = 1 61 ch for the model with the smallest maximum traction of 
1/2 cr . . 

The cohesive zone opening displacements along the cohesive zone length for models with 
different maximum tractions are plotted in Figure 8. There are three cohesive zone opening 
displacement curves associated with three cohesive laws with maximum tractions of 1 / 2cr, , cr. , 

and2cr. , respectively. These cohesive zone opening displacement curves along the cohesive zone 

length are not linear (the shape of the cohesive zone is a cusp). Figure 8 also shows that the 
smaller the maximum traction results in the longer the cohesive zone length. These opening 
displacements are used in Eq. (13) to compute the remote applied stresses. 

The computed remote stresses as a function of the maximum traction of all five CZM models are 
plotted in Figure 9. There are two curves shown on Figure 9 representing two different crack 
lengths. The bottom curve shows that the remote applied stresses for the longer crack length of 
a = 10/^ reaches a constant value for maximum traction greater than one cr,, but the top curve 

for the shorter crack length of a = 1 l ch shows that the remote applied stresses has not reached a 
constant valve for the maximum traction at 4 a c . 

The computed remote applied stresses are then used in Eq. (15) to compute LEFM energy release 
rates G(c r K ). The ratio of G(cr x ) to <D c is defined as the normalized LEFM energy release rate 

(NLERR). If the NLERR is near unity, then the cohesive zone model is equivalent to an LEFM 
analysis, and both models would predict the same failure load. The NLERR as a function of 
maximum traction is plotted in Figure 10. The figure shows that for a long crack length, 
a = 1 01 ch , the NLERR approaches unity for maximum tractions greater than 2 a c while for a 

short crack length, a = l ch , the NLERR reaches unity at a much slower rate (a cohesive law with 
a maximum traction greater than 4 cr, needs be used). 

The NLERR requires high tractions to reach unity for the shorter crack length a = l ch because the 

cohesive zone length relative to the crack length is not small (<sc 1). This large cohesive zone 
length, similar to a large scale yielding zone, may significantly alter the original K dominant 
stress field around the crack tip [13], Hence, the results from the CZM model cannot be 
equivalent to an LEFM analysis. To confirm this hypothesis, the NLERR as a function of the 
ratio of cohesive zone length to crack length is plotted in Figure 1 1 for both the crack lengths. 
This plot reveals that the two curves shown in Figure 10 can be collapsed to become a single 
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curve. This indicates that the ratio of cohesive zone length to the crack length may be an 
important parameter for fracture predictions using CZM. The ratio of cohesive zone length to 
crack length, p / a , must be close to zero for the CZM to be equivalent to an LEFM analysis. 
Physically, this means that regions ahead of the crack tip that exibit mechanisms other than 
brittle fracture (plasticity, bridging fibers, etc) must be very small relative to the crack length for 
LEFM to be valid, and hence, for the CZM approach to be equivalent to LEFM. 

Summary 

A MATLAB® M-file was implemented for numerically solving the CZM integral equations with 
an iterative procedure. Linear softening cohesive laws with a constant CWR but different 
maximum tractions were used for modeling the traction-separation behavior in the cohesive 
zone. Results show that a reduced cohesive zone length is predicted for a cohesive model with a 
large maximum traction, and the ratio of cohesive zone length to the crack size needs be very 
small for the prediction from the CZM to be equivalent to an LEFM analysis. 

The conditions required for a cohesive zone model and an LEFM analysis to be equivalent were 
identified from this study and other published literature. These conditions are as follows: 

1. Stress intensity factor at a fully developed cohesive zone tip needs to vanish if the J- 
integral value around the cohesive zone is set equal to the cohesive work rate. 

2. If the J-integral value around a fully developed cohesive zone is the only parameter to 
determine failure, then the maximum traction of a cohesive law can be changed while 
keeping the CWR constant without affecting the failure load prediction. 

3. The ratio of the cohesive zone length to the crack length needs to be small, so the 
existence of the cohesive zone cannot significantly change the stress field near the crack 
tip. 

4. The maximum traction of a CZM needs to be high. High maximum tractions result in a 
short cohesive zone length. 
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Figure 1 Fracture analysis of a cracked infinite plate using a cohesive zone model. 
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Figure 2 Cohesive zone fully developed at crack growth initiation and unchanged during growth 
[7], 
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Figure 3 J-integral paths around the cohesive zone. 
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Model A Model B Model C 

Figure 4 Five cohesive laws with the same cohesive work rate. 
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Figure 5 Solution procedure of the iterative method. 
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Figure 6 Cohesive zone lengths as a function of maximum traction of different cohesive laws. 
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Figure 7 Cohesive zone lengths vs. crack length for models with different maximum tractions. 
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Figure 8 Cohesive opening displacements along the cohesive zone length for models with different 
maximum tractions ( a = ll ch ). 
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Figure 9 Remote applied stress for models with different maximum tractions. 



Maximum Traction /< t c 

Figure 10 LEFM energy release rates for cohesive laws with different maximum tractions. 
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Figure 1 1 LEFM energy release rates at the moment of growth initiation as a function of the ratio of 
cohesive zone length to crack length. 
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